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Abstract. For subordinators Xt with positive drift we extend recent results on the structure 
of the potential measures 



Jo 



and the renewal densities u^'^^ . Applying Fourier analysis a new representation of the potential 
densities is derived from which we deduce asymptotic results and show how the atoms of the 
Levy measure 11 translate into points of (non)smoothncss of u^'^\ 



1. Introduction 

Let X = {Xt)t>o be a Levy process, i.e. a real-valued stochastic process with stationary and 
independent increments which possesses almost surely right-continuous sample paths and starts 
from zero. If X has almost surely non-decreasing paths it is called a subordinator and can be 
thought of as an extension of the renewal process which accommodates for infinitely many jumps 
on a finite time horizon. In fact any such process can be written as 

= (5< + ^A„ 

s<t 

where A. are random, positive jumps drawn from a Poisson random measure. The non-negative 
coefficient 5 is referred to as drift of the subordinator. The jumps A. are described by a deter- 
ministic intensity measure n((ix), usually referred to as Levy measure, that is concentrated on the 
non-negative reals and satisfies the integrability condition 

/•CO 

(1.1) / {l^x)W{dx) <oo. 

Jo 

The main object of our study are the g-potential measures, q > 0, 

(1.2) U^'^\dx)^E 

of X whose distribution function will be denoted by (x) = J/^'^-' ([0, x]). The slight abuse of 
notation will always become clear from the context. Being the central object of the potential 
theory for subordinators, i/^'^-' is a well-studied object. For more information we refer the reader 
to Chapter 3 of |B96) . Whenever S > each point x > is visited by X with positive probability 
and it is known that in this case the measure admits a bounded and continuous density 
u'^i) = (U'-'^^y with respect to the Lebesgue measure. 

Our main objective is to study the regularity of the g-potential densities u*-') in more depth. 
Under mild restrictions on the Levy measure the main results concern differentiability properties 
of u^"?) when 11 has atoms that may accumulate at zero. We characterize the set of values where 
differentiability of u'-*^ fails as well as the order of derivatives where smoothness breaks down for 
each point of this set. More precisely, if G denotes the set of atoms of Il{dx) then the nth derivative 
of u*-') exists precisely outside of G„ , where Gn contains all numbers that can be represented as a 
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sum of at most n elements of G. Also we provide asymptotic results for u*^') {x) as x goes to zero 
and infinity. For the results at infinity we need to restrict to renewal densities for which from now 
on we use the conventions u = u'^^^ and U = U^^\ 

The paper is organized as follows: Motivation and further background of our research is pre- 
sented in Section [2] An example for the simplest Levy measure 11 = (5i is included for which direct 
computations can be performed to present our abstract results in a simple setting. In Section [3] 
the main results are presented and are proved in Section |4j 

2. Motivation and Background 

2.1. Background. The g-potential measures are the fundamental quantities in potential 
theory and underly the study of Markov Processes in general and Levy processes in particular. 
As subordinators appear naturally in the study of Markov processes (see Chapter IV in |B96) ) 
and are in some sense the simplest Levy processes, it is not surprising that important information 
about its potential measures has already appeared in the literature. For our special case when we 
study a subordinator with positive drift, it is known that with = inf{t > : Xt > x} 

(2.1) u{x)=^-P[Xt^=x], 

i.e. u{x) has the clear probabilistic interpretation as creeping probability of X at level x. From 
this it follows from a renewal argument that u solves the renewal type equation 

6u{x) + / u{x — y)Il{y) dy = 1 







which will be the starting point for our analysis. This simple renewal structure motivates the 
alternative naming of u as renewal density. Here and in the following 



/•oo 

U{y) = U{[y,^)) = / U{dx) 



denotes the tail of the Levy measure. In fact, similar renewal type equations for u^*-' are derived 
below allowing for a study for all q > 0. 

The functions m^''-' are rarely known explicitly but can be described using the universal repre- 
sentation of the Laplace transform of U^'^\dx). On page 74 of B96^ it is shown that 

(2.2) / e-^^C/(dx) = — — , A > 0, 



^(A)^ 

where V'(A) is the the Laplace exponent 

ip{X)^5X+ (1 - e~^^)n(dx). 
^0 

In few exceptional cases such as stable subordinators without drift, i.e. Il{dx) = cx~^~°'dx for 
some a € (1,2), this formula sufficiently simplifies and allows for an explicit Laplace inversion. 

Smoothness issues for the particular case g = have attracted some attention in special cases 
motivated by applications; we shall only highlight two recent results: 

First, it is known that u is completely monotone and henceforth infinitely differentiable for the 
class of so-called special subordinators whose Laplace exponent is a special Bernstein function 
(see |SV06| ). This, as we will show below, is not true for general subordinators. Furthermore, 
for special subordinators, m is a decreasing function. Monotonicity fails in general but we show 
that a similar structural property holds: For all g > 0, is of bounded variation and hence a 
difference of two increasing functions. 

Secondly, when the Levy measure Il{dx) has no atoms and (5 > it has been shown in [CKSlOj 
that u is continuously differentiable. This result can also be recovered by Kingmann's study of 



p-functions (see Chapter 3 of |K72] ) or by our representation (3.4) below. 

We note that the renewal density u is in fact the Kingmann's p-function which is naturally 
associated to a regenerative event, see [K69j and Chapter 3 in |K72j . Using this link one could 
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recover some of our results below for the special case q — 0. However these results are preliminary 
in our study and we focus on finer properties of the functions u'-'?-' for q > 0. 

A further topic of interest are the asymptotics of u and its derivatives. The asymptotic at zero 
and infinity are classically deduced from the renewal theorem. The relation 

(2.3) u{x)^\, as a; -> 0, 



is stated in Theorem 5 in |B96) . whereas 

(2.4) ^''^^^ElXil' 

seems to have been obtained first in |HS01j . We apply two representations developed in the present 
paper to understand the asymptotic behaviour of u'-''^ and (w^''^)' in more detail. It turns out that 
different representations are of very different use: A series representation can be used to deduce the 
asymptotic properties at zero, whereas a Fourier inversion representation is useful to understand 
the more delicate asymptotic at infinity. The use of the Riemann-Lebesque theorem for the Fourier 
inversion representation forces us to assume g = so that those results are restricted to renewal 
densities. 



2.2. Applications. The motivation for our work is mostly theoretical but we outline below some 
applications. Our theoretical motivation stems from the mere importance of potential measures 
in the study of Levy processes. As such those appear in various one-sided and double-sided exit 
problems and especially the probability of hitting a boundary point of an interval at first exit out 
of it is represented in terms of u or if the Levy process is killed at an independent exponential 
time in terms of m^"^^ (see Theorem 19 of Chapter VI in |B96| ). 

Potential densities appear equally in other theoretical studies; For example, representation 



(3.4) below was used to guess the necessary and sufficient analytic condition for existence of right- 
inverses of Levy processes which then, in |DS10j . have been proved via subsequent discoveries. 

Our results should be applicable for simulations of subordinators as detailed knowledge of 
the set of non-smoothness of w'-'^-' could significantly speed up the numerics behind the simu- 
lation/computations of potential densities. This has been noticed in relation with scale func- 
tions by Kuznetsov, Kyprianou and Rivero (private communication) caused by the celebrated 
identity W'{x) = u(x) between scale functions of a spectrally negative Levy process and the 
down-going ladder height subordinator (for more information see Chapter 7 of |B96j ). The 
same holds for potential measures of subordinators. Recall that the scale function is posi- 
tive and increasing, and is involved in computing double exit probabilities for spectrally neg- 
ative Levy processes X, i.e. processes that can jump only downwards. More precisely, with 
T[_a,b] = inf{t >0: Xt>b or Xt < -a} 

^(^T[-y,.] >x)^ w^T^' all 2: > and y > 0. 

Finally, let us point out a possible extension. For atomic Levy measures our main results. 
Theorem [2] and Theorem [sj link the atoms of 11 to points of non-differentiability of w^'^^ . It 
seems possible to extend this to higher order singularities of Il{x). Such a generalization would 
be extremely useful for numerical computations of the scale function which is a basic quantity 
in insurance mathematics as it allows for estimation and calculation the double exit probabilities 
which in the actuarial mathematics context are ruin probabilities. For more background see }CY05) 
and |AKP04j . 



2.3. Example. To motivate our findings, here is an illustrative example for which the renewal 
density u can be calculated explicitly and some non-trivial connections between the atom of 11 
and differentiability of u can be observed. 
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Figure 1. u (solid curve), u' (dashed curve) ioi 5 — 1 and 11 = 5i 



Example 1. Suppose 5 = 1 and 11 consists oj just one atom at 1 with mass 1, i.e. H = Si, 
implying that JI = l[o.i]- For x <1 one immediately obtains 

u{x) = V[Xt^ ^x\= V[no jump before x] = e^^^^") = e"'' . 

For X € [n,n + 1) we exploit the structure of the process more carefully. Since x < n + 1, on the 
event of creeping at x there can be at most n jumps of height one before exceeding x. If precisely 
i < n jumps occur, they need to appear before x ~ i as otherwise the drift would have pushed X 
above level x before the final jump implying 

n 

V[Xt^ = x] = P[exactly i jumps before x ^ i]. 

i=0 

As the number of jumps before x — i is Poissonian with parameter x — i we obtain 
u(x) 



e ^ : a; < 1, 

Er=i^e-(--) :xG[n,n + l) 



The graphs of u and it's derivative are plotted in Figure^^ This first representation already sheds 
some light on the possible behavior of u for atomic Levy measures: in general u is not monotone, 
is asymptotically linear at zero, and possibly non-differentiable. 

Let us consider the issue of differentiability in more detail. Apparently u is infinitely differentiable 
in the open intervals {n,n + 1), n = 0, 1, with 



u'{x) = 



~e : X <l, 



The only problematic points are the integers at which the piecewise defined function u' gets the 
additional summands 

L (n-1)! n\ ^ ■ "■ ^ 

when going from x < n to x > n. This sheds some light on differentiability as for x = I the right 
derivative is 1 and the left derivative whereas for n > 1 the remainder summand vanishes. In 
particular, this implies that u is differentiable everywhere except at 1 with 

u' {x+) - u' [x-) = Ii{{x}). 

Taking higher order derivatives, the same reasoning shows that at {n}, u is {n — l)-times but not 
n-times differentiable. 
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The example reveals three properties which we aim to understand for general subordinators; 
higher order differentiability of u depends crucially on the atoms of H, u is generally not monotone 
and u' vanishes at infinity and is asymptotically linear at zero. 



3. Results 

Our Fourier analytic approach forces us to control the small jumps. To this end we recall the 
Blumenthal-Getoor index /3(n) for subordinators: 



/3(n) = inf |7 > : ^ xm{dx) < ooj. 



Some of our main results arc formulated under the following assumption. 
Assumption (A). Assume that /3(n) < 1. 



Note that property ( 1.1 1 implies that /3(n) < 1 so that the assumption only rules out boundary 
cases and in particular any stable subordinator with drift is included. More direct computations 
for particular examples of interest show that our results will be still valid for Blumenthal-Getoor 
index 1. 

3.1. Series and Integral Representations. The main objective are the potential densities. 
Our results are motivated by the following extension of the renewal type equation for u. 

Lemma 1. The q-potential measure U^'^\dx) has a density u*^"'-' satisfying the renewal type equa- 
tion 

(3.1) 5u'^'>\x) = l- f u^'''>{x-y){Il{y) + q)dy. 



Iterating the renewal type equation (3.1), one heuristically obtains a series representation for 
the potential density u^'^\ Making this a rigorous statement is slightly involved as fl might have 
a singularity at zero. The problems caused by the singularity can be circumvented taking into 



account only the local integrability of 11 at zero which is a direct consequence of the property ( 1.1 ). 
In |CKS10j a proof for the following proposition was carried out for non-atomic Levy measures. 
In fact this proposition can be recovered from Chapter 3 in |K72j where it is disguised in terms of 
p-functions. For completeness we sketch a proof below. 

Proposition 1 (Series Representation for u'--'^^). The series representation 

(3.2) .(^)(.) = 5: y^(i*(ii +,)*")(.) 

n=0 

holds for the potential densities . Here, 1 denotes the Heavyside function. 



At this point we should note that for (7 = 0, Equation (3.2) appears in a slightly different form 
in |UKS10| . This is due to a different definition of convolutions. In this paper we define the 
convolution of two real functions / and g by 

f*g{x)=j f{x - y)g{y)dy, 
Jo 

(and 1 * f*^{x) = 1) whereas in [CKSlOj the convolution was defined by / * g{x) = Jq f{x — 
y)g'{y) dy. In any case, with F{x) = f{t)dt and G{x) = g{t)dt, we have F*G — l*f*g, 
where the second convolution is in the sense of the present paper. 



Unfortunately, deriving properties of u'^'^^ from (3.2 1 is a delicate matter as one has to deal with 



an infinite al tern ating sum of iterated convolutions. Nevertheless, one particular property that 



follows from (3.2) is that is of bounded variation. 
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Corollary 1. There are increasing functions u"-^^ and u^2^ such that 



In particular, it'-'^-' is a function of bounded variation. 

To obtain a deeper understanding, we introduce a more carefull representation for u^'^^ based 
on Fourier analysis. In the following we denote by 



/oo 
e-'-f{x)dx 
-oo 



the bilateral Laplace transform for complex numbers s — X + i9. As n(x) is only defined for a; > 0, 
we put n(a;) for x < 0. 

Proposition 2 (Laplace Inversion Representation of u^'^). Under Assumption (A), for arbitrary 
integer N larger than and any fixed A > 0, 



(A + i9)S^+^{l + i£(n + q){X + 10)) 
is a well defined absolutely integrable function in 9 and 

(3.3) u(')(a:) = * (n + lY^i^) + ^''^ / e^'^s'^'nA + i9) d9. 



ra=0 



There is a good reason to not only consider the simplest case TV = 1. Larger N makes the 
integrand decay faster which is needed to derive higher order differentiability via interchanging 
differentiation and integration. For > 1 it is still easy to work with representation (3.3) as 
the finite sum can be tackled termwise and for the integral classical convergence results can be 
applied. 

As a first application of Proposition [2] we derive a representation for the first derivative of the 
potential densities. In contrast to Proposition |2] we are not free to choose A^ = 1; the minimal A^ 
now depends on the Blumcnthal-Getoor index. 

Corollary 2. The right and left derivatives of u'''^^ exist and are given by 

(3-4) {u^^^)\x+) = + E ^(n + 

(3-5) {u^^^n^-) - -fc^ + E ^(n + 

ri=2 

Moreover, under Assumption (A) the following representations hold for N > ^^j^ijyirx and X > 0: 
(3.6) {u(^^)\x+) = - "("+^ ^ ' + E + + ^'^^ r e^'^h^'^iX + ^9)d9, 



p ^ ^n+i \ • -^1 y ' • - 27r 

n=2 

f -l\n 1 /-oo 



(3.7) {u^'^^Yix-) = - j + + E + + ^'^^ fJ'^h^-{X + ^9)d9 

For every A > the function 

h-^HX + ^0)^ (-An + .)(A + ..))- 



SN+l(l+^C{IL + q){X + i9)) 

is absolutely integrable in 9. 
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Remark 1. We should note that a refinement of the approach of [CKSlOj for non-atomic Levy 



measures yields the series representation {3.4), [3.5) without assuming Blumenthal-Getoor index 



smaller than 1. Indeed, the difference is that in our setting H can have jumps and therefore is 
outside of the scope of |CKS10| . However, the jumps do not affect the study of the uniform and 



absolute convergence of the series in (3.4) and (3.5), and it can be carried out as in |CKS10| 



The arguments of |CKS10| are not repeated here; we only prove the more usefull Laplace inversion 
representation under Assumption (A) and derive from this the series representation. 

Remark 2. We do not know how to derive the asymptotics of u' at infinity only from the series 



representation (3.4-), (3.5). It seems that a representation of the type (3.6), (3.7| is much more 
useful for studying the properties of u' and therefore it is an interesting problem to get such a 
representation without Assumption (A). 

3.2. Higher Order (Non)DifFerentiability for the Potential Densities. Having established 
that u^'^) is of bounded variation in Corollary [I] u^'^) must be difFerentiable away from a Lebesgue 
null set. The null set can be identified as the set of atoms as a consequence of either of the 
representations given in Corollary [2j 

Corollary 3. The potential densities u^'^) are differentiable at x if and only if x is not an atom 
o/n. More precisely, 

(3.8) = 

It is interesting to see that the derivative of u'-^-' only jumps upwards and how the size of the 
jumps is determined by the weight of the atoms and the drift. The reader might want to compare 
this with Figure [T] 

Unless explicitly mentioned, from now on we assume that Assumption (A) holds. 

The remaining part of this section deals with higher order differentiability properties of the q- 
potentials. The study focuses on subordinators whose Levy measure has atoms which accumulate 
at most at zero. To emphasise the effect of atomic Levy measure on differentiability, the results 
are divided into three theorems for non-atomic Levy measure, purely atomic Levy measure, and 
Levy measure with atomic and absolutely continuous parts. 

The case of non-atomic Levy measure already appeared in |CKS10j . For Blumenthal-Getoor 
index smaller than 1 their results follow easily from the Laplace inversion representation. To 
present a complete picture we state and reprove the following result of jCKSlO] for smooth Levy 
measures. 

Theorem 1. Ifti is everywhere infinitely differentiable, then C/^'^ is everywhere infinitely differ- 
entiable. 

More accurate effects of the atoms on differentiability are revealed in what follows. The explicit 
calculation of Example [T] for \i = 5i indeed suggests more interesting behavior for higher order 
derivatives. In that particular case, u is (n — l)-times differentiable but not n-times differentiable 
at n. The critical points n G N in the example have the property that they can be reached by 
precisely n jumps of the size of the atom 1. 

Denote by G the atoms of the Levy measure 11. We say that a; > can be reached by n atomic 
jumps li X — X)"=i 9i with gi € G. For fc > 1 we define the sets 



Gfc = J x > : .X = e G, 1 < J < fc 



of points that can be reached by at most k atomic jumps. The next theorem shows how higher 
order differentiability is connected to the sets Gk- As it is formulated for purely atomic Levy 
measures it can be seen as the counterpart of Theorem [T] 
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Theorem 2. If II is purely atomic with a possible accumulation point of atoms only at zero, then, 
for k > 1, 

is {k + l)-times differentiable at x <^=> x ^ Gk- 
In particular, if x cannot he reached by atomic jumps only, then U'^''^ is infinitely differentiable at 

X. 

The foUowing example shows that we can easily find examples with exotic difi^erentiability 
properties. 

Example 2. // 11 has atoms onl/k for k > 1, then U^'^^ is infinitely differentiable at x if and 
only if x ^Q. For x = n/k, U'^'^^ is at most n-times differentiable at x. 

Finally, we show how Theorems [I] and [2] combine each other: the atoms prevent ?7^'^ from being 
twice differentiable and an additional absolutely continuous part does not change the behavior. 

Theorem 3. If Ii{x) — ni(a:) + I\.2{x), where Hi is purely atomic with possible accumulation 
point of atoms only at zero and 112 is infinitely differentiable, then, for k> \, 

U^'^'^ is {k + l)-times differentiable at x <S=> x ^ Gk- 
In particular, if x cannot be reached by atomic jumps only, then C/*^'?) is infinitely differentiable at 

X. 

Remark 3. In fact Theorem is still true if 1^2 possesses only finitely many derivatives. Then 
the statement will be valid, for any k, such that H2 is k-times differentiable. 

Remark 4. Recalling the close connection of u and creeping probabilities given in \2. 1^ , it would 
be interesting to have a good probabilistic interpretation of non-existence of derivatives of certain 
order for the creeping probabilities - 

3.3. Asymptotic Properties of Potential Densities. Asymptotic properties at zero and in- 
finity of U are classical results in potential theory of subordinators (see for instance Proposition 
1 of Chapter 3 in |B96) ). Also the limiting behavior of the renewal density u at zero and infinity 
are known (see (2.3) and (2.4)). 

In what follows we aim at giving more refined convergence properties based on the repre- 
sentations for u^"^^ and (m'-'^^)'. Some results will be only valid for q = 0. Utilizing the series 
representation of Section |3.1[ we first strengthen the asymptotic at zero. In the following f ^ g 
denotes strong asymptotic equivalence, i.e. lim ^ = 1. 

Theorem 4. For general subordinator with positive drift and n > 0, the following strong asymp- 
totic equivalence holds: 



5]^l*(g + n)*'=(x)- J«)(:.) 

k=0 



^-^l*{q + nr^"+^\x), asx^O- 



The first order asymptotic of (2.3) appears in the theorem for n — and reveals more precise 
qualitative information: 

Corollary 4. The potential density is asymptotically linear with slope — (n(IR) at zero 
iff the Levy measure is finite- If the Levy measure is infinite, then u^'?) decays faster than linearly- 

A simple class of examples are stable subordinators with drift. 

Example 3. // the Levy measure has stable law with a G (0, 1), the strong asymptotics at zero 
are given by 

1 ^ ^ c 



S Jo y'^ d^l-a)' 

Levy measure putting mass precisely to one point provides an example of asymptotically linear 
behavior. 
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Example 4. // the Levy measure has only one atom at 1 and 6 — 1, then the direct calculation 
in Example^ or Corollary^ both lead to the linear asymptotics 

1 — u{x) = 1 — e~^ X, as X 

As for absolutely continuous Levy measures the existence of the derivative of the renewal density 
was established quite recently, the asymptotics of u' seem to be unknown. On first view one might 
be tempted to differentiate the asymptotics of u to derive the asymptotics of u' . Due to lack of 
knowledge on ultimate monotonicity of u and u' we cannot apply the monotone density theorem 
and instead work with the representations of Section |3.1[ 



Theorem 5. Assume condition (A), i.e. /3(n) < 1. Let 11(0+) = oo, then there exists n > 1 
such that /3(n) < and then 

(«(^)) {x+) = + ... + + 0(1), as X 0, 

t (aWr ^ n(a;-) (-l)"n*"(a;) 

(^W)(^_) = + ^ ^ +0(1), asx^Q. 

//n(0+) < oo then, as x -> 0, 

{x+) = ^ ^ ^+o(l), 

(x-) = ^+o(l). 

For a large class of processes this implies that the asymptotic of {u'^'^'^ Y are indeed given by the 
derivative of the asymptotic of m^'^-' : 

Corollary 5. // 11 has Blumenthal-Getoor index /3 < 1/2, then (u'^^')'(a;+) ^ — ^^^'^fctl (jjjjf 
(u(?))'(x-) - _(9±nK£^ a< zero. 

The hypothesis of the theorem is not necessary to have the asymptotic Il{x)/5'^ at zero as the 
next example shows. 

Example 5. If X is stable with index a € (0, 1) and drift 5 > 0, then 

n*^x) ^Ci f {x- dy = Cs-t-^^+i 

"'0 

and inductively n*"(x) — C„a;~""'^"~^ . This, combined with Theorem^ implies that 

, n(x) Cix-" 

u ix) ~ -5— = -5 — , as a: ^ 0. 



Finally, we consider the asymptotics of u' at infinity which is technically more involved as we 
could not derive the asymptotics from the series representation (3.4), (3.5). Instead, the proof is 
based on a refinement of the Laplace inversion representation. 

Theorem 6. //E[Xi] < 00, then 

lim u' {x+) = lim u' {x—) = 0. 

As our approach does not work for infinite mean subordinators, it remains open whether u' 
vanishes at infinity or not in this setting. 



4. Proofs 

4.1. Series and Integral Representations. 
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Proof of Lemma [7| When q — the statement was aheady mentioned in Section |2.1| and is due 
to Kesten (see (K69| ). Let us now fix g > 0. First note that for any random variable Cg ^ Exp{q) 
which is independent of X using the Markov property at time Cq and integrating out the possible 
positions of Xg^ < x we get 

U{x) = e( /" ' l{x,<.}di) + / U(x~ y)P(Xe, G dy) = U<^'^\x) + qU * u<^'i\x), 



since m^'^-' is continuous and bounded whenever the subordinator possesses a positive drift, see (iii) 
in Theorem 16 in [B96]. Therefore, differentiation yields 

u{x) = u'^'^^x) + qu* u^'^''{x). 

Comparing the latest with 

Su{x) = P(Xt, = a;) 

= V{Xt^ ^x;T,< eg) + f P(Xt(,_„) - x - y)nXe, G dy) 



= P{Xt^ ^x;T,< eg) + q6 / u{x ~ y)u^^'> {y)dy, 

Jo 

we get the expected relation 

V{XT^=x;T.,<eg) = Su('^\x). 
We use this final relation to write 

P(T, < e,) = 1 - P(Xe, < x) = + P(Xt, > x; < e,). 

Next we have, using the compensation formula from Chapter in jB96| . 

V{Xt^ > x; < eg) = q f V l{x_<.. x,>.}' 

Jo ^ ^ 



s<t 

/•OO 



I / e-"' / U{x- y)V{X, G dy)ds 
Jo Jo Jo 

Ilix-y)u<^'^Hy)dy. 

This, combined with P(Xe^ < x) ^ q u^*) {y)dy, proves the assertion. □ 

Proof of Proposition^' Since 11 may not be continuous, we cannot directly apply the general 
Theorem 6 of ICKSlOj . Nonetheless, the proof of Theorem 6 of |CKS10| can be adapted to this 
situation by noting that in fact it is not crucial that g there is continuously differentiable but the 
fact that it is almost everywhere differentiable with a negative derivative vanishing at infinity. The 
latter is due to the fact that g' is used in convolutions. In this case, for completeness, we sketch 
a proof but we refer to [CKS10| for details. First note that each summand of 



(4.1) 0(x) 



n=0 



is finite as for subordinators xll{dx) is finite showing that 11 is not exploding too fast at zero. 
We now show that the series is absolutely converging in particular showing that it is well defined 
and can be rearranged. First, since 1 * (ll + q)*^" ^\x) is increasing we get iteratively 

(4.2) 1* (n + (7)*"(a;) < (l * (ll + g) (:r))". 

As 1 * (n + q){x) is continuously increasing and 1 * (II + g)(0) =0 there is 6 > such that for all 
X <b, 1* {n + q){x) < 6/2 and hence the series defining </) is absolutely and uniformly convergent. 
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Now we use an iterative procedure to extend the absolute convergence to all a; > 0: reasoning as 
above we obtain 

l*(n + g)*"(2&) _ /■''l*(n + g)*("-'^(26-y)(n + g)(y)^^ 



l*(n + g)*^""''(2fe-y) (n + g)(y) 



dy 



b 



(5" S 



^ l*(n + g) ^ '{2b)l*{n + q){b) l*(n + g) ^ '{b)l*{U + q){2b) 
~ (5" S ^ (5" 6 

<(n + l) A*(n + g)W y"''^ l.(n + g)(2b) ^ (n + l)l.(n + g)(26) ^ 
~ \ 5 / i5 ~ 2"~i (5 

Iterating this arguments as in the proof of Theorem 6 of |CKS10| shows locally uniform and 
absolute convergence of the series for all x > 0. To verify that ^ equals u we exploit Fubini's 
theorem to obtain 



n=l 



implying that is a solution of Equation (3.1). Uniqueness of solutions for locally bounded 
functions follows since for two solutions / and cf) 

\f~C^\ix) < \{f-q^)*l/S{n + q){x)\ 

-|(/-0)*l/<5^(n + g)*'(:r)| 
<sup|/-0|(s)l/<5'=l*(n + g)*'=(a;) 

for all k. The right hand side converges to zero since ^"'^l * + qY*'' {x) goes to zero for all a; > 
as shown above since it is a member of uniformly convergent series. □ 

Proof of Corollary [7J The proof follows directly from Proposition [l] To define the functions u^^'^ 
and u^'^ we separate the positive and negative parts as 



1 /i=t , \*2n/ N oo 1 /-pr , \*(2n+l), , 
n=0 n=0 



which is possible taking into account the absolute convergence of the series representation. Each 
summand of the defining series for u^'-* and u'"2 ^ is increasing, thus, u^-^'^ and u'^'^ themselves are 
increasing. □ 

We are now going to derive the Laplace transform representation of Proposition [2] As men- 
tioned above, the main difficulty is that we need to deal with convolutions of possibly unbounded 
functions. Let us first motivate the approach. Taking into account the series representation we 
divide u^'^' into two parts: 

n=0 n>N 

-E^i*(n + '^)*» + ^"''(-)- 

The goal of the following Fourier analysis is to derive a convenient integral representation for 0^'''. 
For this sake first note that it follows as in the proof of Proposition[l]that (j)^''^ is the unique locally 
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bounded solution of the following renewal type equation: 

(__T\N - AT ^ 

(4.3) cj^'^^^x) ^ * (n + ~ * (n + q){x). 

If we were allowed to turn convolution into multiplication in Laplace domain for Re{s) > 0, 
Equation (4.3 1 transforms into 



(4.4) 



1 



/:0^''(s)£(n + g)(s) 



Solving (|44| for Ccj)^^'', leads to 
(4.5) 



whenever the quotient is well-defined. If furthermore we were able to verify integrability conditions 
needed for Laplace inversion we obtain an integral representation for 0^''^. We are now going to 
check what is needed to turn this formal approach into rigorous statements. 

For the rest of this section we set 5 = 1 in order to simplify the notation. 



For the first step it is shown that indeed Equation (4.3) turns into Equation (4.4). A priori this 
is not clear due to the possible singularity of 11 at zero. 



Lemma 2. There is Xo > such that dp^'"^ solves (4-4^ on {X^ iO : \ > Ao, e M}. 

Proof. To show that the first transformation can be carried out we show for s — X + iO, X bounded 
from below by some Ao, that C{n + q){s) and £0^'«(s) are well defined to deduce £(n + g)*"(s) = 
/:(n + g)(s)", £(n + g)*<^^'9(s) =£(n + q)(s)/:0^'9(s), and£l*(n + g)(s) =£l(s)£(n + g)(s). 

To validate Laplace transformation of (ft + g)" for A large enough, note that we may choose Aq 
such that 



(4.6) 



^°''{n + q){x)dx < 1 



which is possible as Jp n(a;) dx < oo and lima;_j.oo n(a;) = 0. It now follows directly from Fubini's 
theorem that iterated convolutions of n + g turn into multiplication under Laplace transforms. 
We now show that ip^'"^ can be Laplace transformed for which we use the fact that 

oo oo 
n=N 7i=N 

to justify the change of summation and integration in the following: 
(4.7) 

/■OO OO oo nOC 

C(t)^^i{s)^ e-*" ^(-l)"(n + (7)*"(:E)dx= ^(-1)" / e" '^li + q)*"" {x) ds < oo. 



N 



Now it is_only left to show Lcp^-i * {Ii + q)(s) = /;</)^''?(s)/:(n + g)(s) and further CI* {Ii + q){s) = 
Cl{s)C{Ti + q){s). As n + g and 0^''' can be Laplace transformed, we obtain from (4.6 1 and (4.7 1 
the bound 



oo /"OO 



Jo 

oo />oo 



cP''-''i{t-x){n + q){x) 



dtdx 



e-^(*-"V^'9(i-x) 



e-^^\l + q){x] 



dtdx < oo 



enabling us to apply Fubini's theorem once more to obtain £ip'^''^*{Il+q){s) = £(l)^''^{s)C(n+q){s). 
The final identity £1 * (II + q){s) = £l(s)£(II + q){s) follows from similar arguments noting that 
|£1(5)| < i. □ 
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The second step in our analysis consists of showing that the convolution Equation (4.4) can 
indeed be solved in Laplace domain leading to Equation (4.5). The following lemma is stronger 
then the previous as we can show that g^''^{s) is well defined for any s = X + i9 with A > even 
though a priori we do not know that g^'' is the Laplace transform of 0^-''. 

Lemma 3. Suppose that 11 is non-trivial, then g^''^{X + iO) is well-defined for A > 0. 

Proof. To show that g^'''{s) is well-defined at s G C with positive real-part, it suffices to show 
that 

C[n + q)(s)^-l. 

Let us assume that £(n -I- q){so) = —1 for some so = Aq + i6o with Aq > 0. Without loss of 
generality way may assume 6*0 7^ as otherwise the contradiction follows trivally. The assumption 
necessarily implies that 

(4.8) Jm(£(n-f g)(so)) = / 6'^°'' sm{9ox){IL + q){x)dx = 0. 

Jo 

As n is decreasing, we see that 

(4.9) / e"-^"'=(n + q){x) dx < 00. 



Dividing the integral of the absolutely integrable function e -^"^(n + q){x) sin(6'ox) into pieces of 
the length of one period of the sine function and applying Fubini's theorem, we obtain 

/•oo 

= / e~^°''{U + q){x)sm{9ox)dx 
Jo 

/•oo ^ 

" / ^^[2k7r/0o,2{k+i)7r/eo)ix)e~^°'^(n + q){x)sm{9ox)dx 

00 2(k + l}^ 

= y / '° e-^'"'{fl + q){x)sm{eox)dx. 

k=0'' 

As e^'^°^n(a;) is strictly decreasing, unless (H -I- q){x) ~ and (fl + q)(x) is non-increasing each 
summand 

e-^°="(n-f^g)(a;)sin(6'oa;)da; 

must be non-negative and hence vanish as the total sum is zero. In particular, this implies that 
n(a;) -I- g = for all a; > so that for g > a direct contradiction occurs. For g = the 
contradiction occurs as H was assumed to be non-trivial. Thus .g(s) is well-defined. □ 

Remark 5. If furthermore ^\X{\ < 00, then g^'^ is well-defined also on the imaginary axis. This 
follows from the same proof noting that in this case (4.9) holds as well for Aq = 0. Indeed as 

2(fc + l)7r 

each term Jatx " ^{x) shi(9Qx) dx has to vanish we conclude that H has to be concentrated on 
{2fc7r/0o}*:>i • On the other hand, in this case, as 

/•oo 

Re{£U{so))= e-^'>^ cos{9ox)Il{x)dx 
Jo 

ynC-^) ['^co.i9ox)dx^o, 



k>0 

we see that i?e(£n(so)) 7^ —1 and thus g'^'''{s) is well-defined. 
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The third step of our derivation of an integral representation for is an inversion approach 
for g^''^. We now briefly discuss the connection to Fourier transforms which is crucial for the 
inversion: for integrable functions / define for a; S M 



/oo 
-oo 



Apparently, the Fourier transform T appears when evaluating the Laplace transform on the imag- 
inary line only. Defining the auxiliary function 

the simple connection is 

Taking into account this close connection of Laplace and Fourier transforms, classical Fourier 
inversion for A > gives the inversion formula (also known as Bromwich integral) 



1 

1 1'°° 



27r 

if (?^''(A + iO) is absolutely integrable with respect to 6. To prove the needed integrability we 
start with a simple estimate. 

Lemma 4. For any a > and y < a the estimate n(j/) < C {a, e)y~^^^^'>~^'^^ holds for all s > 
with /3(n) + e < 1. 

Proof. First note that by the definition of the Blumenthal-Getoor index j/'5(n)-i-£ n(dy) < oo for 
any £ > 0. The claim follows from the simple observation that for any a > there is ^ > such 
that for any t < 5 

a>J y^^^^+'IL{dy)>T^^^^+%n{T)-n{S)) 

as ft is decreasing. Letting r go to zero, we deduce limsup^^Q T''(^'+^n(T) < a. □ 

The need for Assumption (A) comes from the following lemma and its consequences. 
Lemma 5. For any A > and e > the following estimate holds: 



|£n(A + i6»)| < 



Ci : \9\ < 1, 



where C = C{e) > 0. 
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Proof. We estimate the imaginary and real part of £11 separately. For the imaginary part we first 
estimate for 6* > 0: 

\lm{CU{X + i9))\ 

sm{y)rx{y/9)dy 







E 



(rxiy/O) - rxiiy + 7r)/0)) sm{y)dy 



2k-n 



< 



< 



TTy2rx{2k/e)~rx{{2k + 2)/9)+e-^ / 
fc=l ^0 



{rx{y/e)~rx{{y + 7r)/d))dy 



rx{y/9)dy 



+e-l 



where we have used Lemma [4] and that r\ is decreasing in the last inequality. Unfortunately, 
this uniform in A upper bound is not suitable for all 9 as the constant of Lemma |4] explodes as 
6 approaches zero. To circumvent this problem we derive a different upper bound that works 
everywhere equally well but is not uniform in A: 

\lm{CIl{X + ie))\ 



< 



r\{y/0)dy 
rxiy/e)dy^ 



r\{y/d)dy 



< 



< 



Il{y / d)dy + d-^n{l) / e-^y^/'^Uy 



c / {y/e)-^f^(^^+'^dy + e-^i{{\)- 



^C + C-, 

where we again used Lemma |4] but now y/9 does not explode for small 6 as we only integrate up 
to 6. Having an estimate for positive 6 we note that Im{£Il{X + id)) as a function of 6 is odd to 
deduce that 



(4.10) 



\im{cn{x + ie))\ < 



1^1 <1, 

\e\>i. 



Similarly, we estimate the real part 
\Re{Ctl{X + i6))\ 

cos{ey)rx{y)dy 



cos{y)rx{ye-^)dy + e 



= l"'(4fe+l)i 



(rxiye ^)-rxi{y + n)e ^))cos{y)dy 



< 



rxiyO ^)dy 

for large \9\ and precisely as above for small \9\. This finishes the proof of the lemma. 
The upper bound can now be used to derive the necessary integrability of g'^''^. 



□ 
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Lemma 6. For arbitrary integer N larger than and any X > we have 



(4.11) 



\g^''^{X + i9)\ d0<oo 



for g^''^ defined in {4-5). 



Proof. As we have already found a good upper bound for CIl{s) in the previous lemma, it suffices 
to show that the denominator 

p{x + i0) = i + c{u + q){x + ie) 

is bounded away from zero. In Lemma |3] we have shown that p(X + i6) has no zeros for A > and, 
hence, by continuity of p it suffices to show that as \0\ tends to infinity p{X + iO) stays bounded 
away from zero. To this end it suffices to note that from Lemma [5] 

q 



lim L{iv + q){x + ie)= lim {ai{x + ie) , ^ 



= 0. 



Using the fact \Cl{X + id)\ = |l/(A + i6 
the upper bound 



< 



min { j^} for A > 0, we employ Lemma jsj to obtain 



g^'?(A + ie)\< c\ci{x + ie)\\c{ii + q){x + ie)\ 



N 



(4.12) 



= C 



X + it 



/:n(A- 



X + i6 



N 



< 



f~ii 1 

C'|6/|A'(/3(n)+£-i)-i 



<1, 
> 1. 



The right hand side is integrable in Q as by assumption /3(n) < 1 and e can be chosen sufficiently 
small so that Af(/?(n) + e - 1) < 0. □ 

We are now in a position to derive the Laplace inversion representations for u^'^ . 

Proof of Proposition^ For A > Aq, Aq satisfying Jp°° e^'^''^(n + q)(x) < 1 we can directly follow 
the strategy explained before Lemma[2] The proof then follows directly from the definition of ^^''^ 
and Laplace inversion justified by Lemmas [2] [3) and [6) 

The proof of the proposition is complete if we can show that for arbitrary < A < Aq 



r(A) 



(s) ds ■ 



r(Ao) 



where r(A) = {X + i9 : 9 & M}. Laplace transforms are analytic (see for instance Theorem 75.2 
of |K88| \ and g^'^ has no singularity for A > by Lemma [sj hence, Cauchy's theorem applied to 
the closed contour formed by the pieces 

r(Ao)n{|e| <i?}, 
T{X)f^{\e\<R), 

$(i?) = {s: s = r + iR,re [A,Ao]}, 

$(i?) = {s: s = r-iR,re[X,Xo]}, 

taken with the right orientation implies the claim. Note that the integrals over the horizontal 
pieces vanish as R tends to infinity: 

e'''g^''^is)ds <Ce^°^(Ao-A) lim |i?|^'(/3(n)+e-i)-i^ 



lim 



_R-i.oo 



where we have used (4.12) for \9\ > 1 to estimate \g'^''^{s)\ for R big enough. Choosing e small 
enough so that /3(n) + e — 1 < 0, the right hand side tends to zero. The same argument shows 
that the integral over $(i?) vanishes. □ 

Proof of Corollary As remarked after the corollary, the arguments of |CKS10) will not be re- 
peated. Instead, under Assumption (A) we prove the Laplace inversion representation and deduce 
from this the series representation of |CKS10| . 

We first show that the right and left derivatives of u^^^ exist and are given by the representation 
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of the theorem. First, right and left derivatives of the finite sum in ( |3.3| exist by termwise dif- 
ferentiating the finite sum and using that 1 * (fl + (7)*"(a:) — (fl + q) "(y) dy is differentiable 
from the left and the right with derivative (fl + qf^ix—) (resp. (fl + g)*"(a;-|-)). As iterated 
convolutions are continuous, only the first summand is not everywhere differentiable. 
To see that the integral is differentiable at x and to deduce the integral representation of (u*^^')' 
note that 



dx 2iT 



e'''''g^'''{X + i9)d9 - 



1 

1 

2^ 



{X + i9)e'^''g^-''{X + i9) d9 
e*^^/i^''?(A + i6') d9. 



The differentiation under the integral is justified by dominated convergence and the upper bound 



dx 



^^g^^'^{\ + i9) 



f^''{\ + i9) \ < 



C\9\j^ 
C|6»|JV(^(n)+e-i) 



1^1 <1, 
\9\>1, 



derived in (4.12) which is integrable in 9 for sufficiently small e by our choice of N. 

As a second step we now deri ve the pointwise series representation from the Laplace transform 
representation of (w^^^)'. As for (4.121 we obtain the upper bound 



(4.13) 



\h'^-'^{X + i9)\ < C\CU{X + i9)\ 



N 



< 



C\9\ 



N{Pin)+e-l) 



<1, 
> 1, 



for arbitrary £ > 0. To prove the corollary it suffices to show that for N tending to infinity, the 
Laplace inversion integral 

1 



2tt 



eiS^h^^il^X + iff) d9 



vanishes for fixed a; > and A > 0. With our choice of e, i.e. /3(n) + s — 1 < 0, (4.131 impfies 
pointwise convergence e*^^/i^''^(A + i9) 0. We are done if we can justify the change of limit and 
integration. This comes from the uniform (for N > 2[(/3(n) + e — 1)^^] + 2 ) in 6* integrable upper 
bound 



'^^h^'i{X + i9) \ < 




\0\ < 1, 
1^1 > 1> 



and the dominated convergence theorem. 



□ 



4.2. Higher Order (Non)Differentiability. In this section the results on differentiability are 
proved. In contrast to Corollary [3) which follows either from differentiating the Laplace inversion 
representation or differentiating termwise the series representation of u^'^\ the proofs for higher 
order derivatives are exclusively based on the more elegant Laplace inversion approach. This 
forces us to assume 13 {IV) < 1 and we do not see how to circumvent this (probably dispensable) 
restriction. 

To reduce the proofs to (non) differentiability of iterated convolutions, differentiability of the 
Laplace inversion integral is ensured in the next lemma if N is sufficiently large. 



Lemma 7. For N > 



jj^— j-, the Laplace inversion integral in [3.3) is everywhere k-times 



continuously differentiable in x with 



1 



da;*-' 27r 

Proof. To check that we can differentiate 



1 1"^ 

^SxgN,q^X + z0) d0 = — / e(^+"')"(A + i9)''g^'''iX + i9) d9. 
2tt I „ 



1 



27r 
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/c-times under the integral, by dominated convergence we need to find an integrable upper bound 
for the derivative: 



dx'^ 



■\k N.qi 



Integrabihty follows directly from the choice of N and the integrable upper bound (4.12) for 
sufficiently small e. □ 

Having understood differentiability of the remainder term, to prove higher order differentiability 
of u*-') we choose N large enough to apply the previous lemma and then deal with the finite sum 
of convolutions. Here is a lemma for smooth Levy measures. 

Lemma 8. // / : M — )• with f{x) — for x < is infinitely differentiahle on and locally 
integrable at zero, then /*" is everywhere infinitely differentiahle and integrable at zero for any 
n>l. 

Proof. The proof is easily conducted by induction with basis n = 1, i.e. f{x). Then the simple 
identity 

f *^^(:c-y)f{y)dy+ f*- {y) f {x - y)dy 



and the induction hypothesis confirm the statement of the lemma with respect to differentiability. 
The integrability follows similarly from the representation 

' f<^+^\x)dx^ C f*'\y) C ' f{s)dsdy 







and the integrability of f{x) at zero. □ 

Combining the previous lemmas we can prove Theorem [l] 

Proof of Theorem^- As the potential measure J7*^'^ is differentiahle with derivative u'^'^\ to differ- 
entiate (fc + l)-times U^'^^ it sufhces to differentiate fc-times the potential density u'^'^\ Applying 
Proposition j2j with N > — ^(n)-i yields 

j/c+l Jfe 



dx^+^ dx^ 
_ d^ 
dx^ 



As N was assumed to be large enough, the integral is everywhere fc-times differentiahle by Lemma 
[7) For /3(H) close to 1 the sum might be arbitrarily large but is always finite. So we can differentiate 
termwise the sum and use ^1 * (H + q)*'^{x) = (ft + q)*''^{x) as H is continuous leading to 

(4.14) 

-,U^'\-)=Y. W;^(n + g)™(.) + e^-- / e'^{\ + ^efg^■'^{\ + ^B)de. 



n—1 



oo 



Applying Lemma [8] to each summand of the finite sum concludes the proof. □ 

The proof revealed the full strength of the Laplace inversion approach compared to the series 
approach of jCKS10| . Their major technical problems consist of justifying differentiation under 
the alternating infinite sum (3.2). As we split the infinite sum into a harmless finite sum and an 
integral which can be dealt with easily, the main problems of [CKSlOj have been circumvented. 

Next we analyze the influence of atoms of H on (non) differentiability for which we start with a 
lemma on higher order convolutions for discrete Levy measures. 

Lemma 9. Suppose H is purely atomic with possible accumulation point of atoms only at zero, 
then 

(a) H*" is a polynomial of order at most n—1 away from Gn for n>\. 
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(b) n*" is everywhere (n — 2) -times dijferentiable but not (ji — 1) -times differentiable at 
Gn\Gn-i forn>2. 

Proof. The proof is based on the simple observation 



(4.15) l[0,a]*/W 



Iq fiy) dy ■■a>x, 



for integrable / vanishing on the negative half-hne. Hence, l[o,a] * / is continuous everywhere and 
differentiable at x and x + a if and only if / is continuous at x. We will resort to the fact that as 11 
is purely atomic, 11 = Tl{x, oo) is piecewise constant and can be represented as linear combination 
of step functions, i.e. 

(4.16) n(.T) = ^n({a})l[o,a](a;). 

aeG 



To prove (a) we proceed by induction which we start with n = 1. Taking into account (4.16), 
n*i = n is a polynomial of order away from G. Next, assume that n*" is a polynomial of order 
at most 71 — 1 away from G„. To prove the claim for n*^""'"^^ we fix an interval {d,b) such that 
d e Gn+i U {0}, b e Gn+i and (d, b) n G„+i = 0. Our aim is to show that for every x € (d, b) there 
is an open interval A(a;) C (d, b) such that n*("+i) is polynomial of order at most n on A(a;). Fix 
X G {d,b) and choose a*{x) = maxjc € G U {0} : x — c > d}. As the atoms accumulate only at 
zero there is a neighbourhood A(a;) of x such that z — a*{x) > d for every z e A(x). Next, by 
Fubini's theorem we observe 

(4.17) n*("+i)(a:) = J2 n({«})i[o,a] * n*"(x). 

aeG 



By the induction hypothesis and (4.15) all summands are locally polynomials of order at most n. 



To show that n*^""'"^) is a polynomial of order at most n, we split H*" according to the two cases 
of ( |4l5l ): 

n*("+i)(:r)== J2 n({«})i[o,a]*n*"(x)+ J2 n({a})i[o,,] *n™(:r) 

= ru*"{y)dy J2 n({a})+ "({a}) f fl*'\y)dy. 

The first summand clearly is a polynomial of order at most n on A{x) by the induction hypothesis 
using that X]a>a; n({a}) is a constant on A(a;). To show that the second summand is a polynomial 
as well, we write for all z S A(x): 

J2 n({a}) r U*-iy)dy 
= J2 ma})£ n*-iy)dy+ n({a})£^n™(y)dy. 

a* {x)>aeG a* {x)<a<z,aeG 

We clearly deduce by the induction hypothesis that the second sum, having finitely many sum- 
mands only, is a polynomial of order at most n on A{x). By the definition of a*{x) and the 
induction hypothesis, 11*" is the same polynomial on {d,b). Also, since for z — a > d 

I lV"'{y)dy<a sup |n*"(s)| < Ga, 

the first sum is clearly absolutely sumniable and henceforth it defines a polynomial of order at 
most n. Thus we conclude that on A{x) we have that n*("+i) is a polynomial of order at most n. 
Representing (6, d) as a union of neighbourhoods A(a;) shows that n*'^"+^^ is indeed a polynomial 
of order at most n on (fe, d). 

In particular, (a) shows that IT*" is infinitely differentiable away from G„. To prove the claimed 
non-differentiability in (b), a different approach is needed. We prove the assertion again by 
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induction in n. The first step is to show that H*^ is everywhere continuous and not difFerentiable 
at G2\G. Continuity foUows from the continuity of l[o,a] *n(a;) and the locaUy uniform convergence 



of (4.17) which can, applying monotonicity of 11 and the properties of H, be seen for e small enough 



and N large enough from 

sup[ n({a})i[o.a] *n(x)+ n({a})i[o,a] *n(a;) ) 

^ VoaSG N<aeG ) 

= sup( Y n(W) r n(2/)dy+ Y n({«}) f nv)dy 

^ \,^„^r. Jx-a N<aeG "'0 ' 



< sup ( Y n({a})an(x ~e)+ f n(y) dy Y ^({a}) 
<supC(x)( Y n({a})a+ ^ n({a}) 



N<aGG 



N<aeG 



E-> 0,jV^ OO _ 



where C(x) are locally bounded constants. 

Now we show that 11*^ is not differentiable at G2\G. Although we already know that n*^ is a 
polynomial away from G2, we derive a second representation showing that for a; ^ G2 it can be 
differentiated termwise. As the atoms are discrete, there is a largest atom c < x implying that 



Gn (c, x] = 0. Taking into account (4.15 1, that 11 is constant in (x — c, x], and the definition of G2, 



we see that l[o,a] *n is infinitely differentiable away from G2. To show that the infinite sum (4.171 
can be differentiated termwise for n = 1, we derive a locally absolutely and uniformly summable 
in x upper bound for the series of derivatives: 



En(W) 

oGG 



^l[0,a]*n(.) 



Y n({a})n(x)+ Y n({a})|n(x)-n(x- 



x<aeG x>a£G 



The first summand is bounded as n is locally constant. For the second, the sum only runs over 
X — c < a as otherwise n(x) — Ii{x — a) = 0. There is no accumulation point of atoms at x, so 
the second summand is bounded by Cn([x — c,x\) < 00 and, hence, fl*^ can be differentiated 
termwise. 

With this in hand we can show non-differentiability at G2\G. Let b g G2\G, then, since G has a 
possible accumulation point only at 0, there is a neighbourhood A{b) of b such that [y, a:]nG2 ~ {b} 
for y,x e A{b). Since b G G2\G we have n(z) = U(b) for all z G A^b) so that 

(n*2)'(6+) _ (n*2)'(5-) 

= lim [(n*^)'(x)-(n*2)'(y)" 

xlb, y1;b L ^ ' ^ ' . 

= JrS ^ n({a})(n(x)-n(y))+ Y n({a})(n(x)-n(x-a)-n(?y) + n(y-a))) 



b<aeG b>aeG 



= Y n{{a})n{{b-a}). 

b>a£G 

Since b G G2\G and G has a possible accumulation point only at zero, the last sum has a finite 
number of summands and is strictly positive. The exchange of limit and summation is possible 
since for all y and x sufficiently close to 6 and all a G G, such that a < c, for some c > 0, 
y — a € A{b) and x — a G A(b). The latter implies that Il{x — a) = n(y — a) and the sum in the 
limit is a finite sum. This proves (b) for n = 2. 
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Next, assume that 11*" is everywhere {n — 2)-times differentiable and that for b e G„\G„_i 

jn-l _ jn-l _ 

< .n*"(6+) - .n*"(6-) < oo. 



To show that at the critical points n*("+^'' is everywhere (n — l)-tinies differentiable we again 
verify termwise the diflterentiability: 



En(W) 



Jn-2 



E n(W) ^(n*"(.)-n-(x-a)) + E n({a}) 



x>aeG 



x<aeG 



Jn-2 



n*"(x) 



As n*" is everywhere {n — 2)-times continuously differentiable, the second sum is locally bounded 
by the property (1.1). Clearly for any x E G„\G„_i there is an interval (x — d,x) such that 
{x — d,x) n G„ = and d < x because G has a possible accumulation point only at zero. The 
latter also implies that there are at most finitely many atoms a such that x > a > d and therefore 
we need to study only 



E n(W) 



d>a£G 



— (n*"(.)-n-(x-a)) 



But from (a), 11*" (y) is a polynomial of order at most n — 1 for all y € {x — d,x) so that by the 
mean value theorem 



(4.18) 



jn-2 



da;"-2 



n*"(x) 



711-2 



dx"-2 



n*"(a;-a) 



< Ca. 



As X]a<d^n({a}) < oo we can interchange differentiation and summation and then using the 
induction hypothesis that H*" is (n — 2)-times continuously differentiable everywhere. In total we 
conclude that n*("+i) is {n - l)-t imes continuously differentiable everywhere. 



Finally our task is to show that for b e G„-|_i\G„ 

d" 



o< _n*("+i)(6^ 

da;" 



dx" 



n*("+i)(6-) < oo. 



As b E Gn+i\Gn and there is a neighbourhood A{b) of b such that n G„+i = {6}, the 
arguments above imply that we can differentiate termwise n-times (4.17) on A{b)\b to get 



f^n _ in _ Jn 

n*("+i)(x)^ E n({a})— i[o,]*n*"(a;)+ E n({a})— 1[„,,] *n*"(a 

a;-a^/l(6),aeG a;-aeA(b),aeG 



da; 



Using (4.15) we rewrite the latter as 



d" 
da;"' 



-n 



*(n+l) 



(a;) 



jn— 1 



= E n(W)^n*"(x)+ E n(W)(;^n*"(a;)-;^n*"(x-a 



a;<aeG 



x>aeG 



jn — 1 



da;" 



da;" 



Now, as 6 € G„+i\G„, n*"(a;) is continuously differentiable on A{b), n({5}) — and on A{b) 
dx^-i n*"(x) = n*"(y), because according to Lemmak)! n*"(a;) is a polynomial of order at 
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most n — 1 on A{b), we obtain 

dx" ^ ' da;" 

= lim [ - — 



in jn 

^ ^ da;" ^ ' 

lim ( ^n*("+i)(a;) - — n*("+i)(y) ) 
,&.yt6 V'^^;" ^ ' da;" / 

lim y 4 -ir'\x -a)- 4 T^*"Hy 

xib.ytb ^ \dx"-^ ^ ' da;"~i 



^ n(W)(|::;?rn™((^-«)+)-£in-((^ 

h>aeG ^ 



Since h G G'„+i\G'„ we have that either h — a ^ Gn or & — a G Gn\Gn-i, where the latter is 
possible only for finitely many, but more than zero, a G G. We mention that the interchange of 
limit and summation is valid since for all a; and y sufficiently close to 6, there is c > 0, such that 
for all a G G such that a < c, x — a G A{b) and y — a £ A{h). We have already argued above that 
the n — 1-st derivative of 11*" is a constant on A{b) due to Lemma [o] By the induction hypothesis 
for those a 

Thus we conclude the induction and the proof of the lemma. □ 

The observation of the lemma motivates the strategy for the main proofs: first, choose N large 
enough that the integral of the Laplace inversion representation can be differentiated as often as 
needed. Secondly, consider the finite sum of iterated convolutions for which critical points for fcth 
derivatives only occur for the fcth summand. Lower order convolutions vanish and higher order 
convolutions are smooth enough. 

Before we give the main proofs, one more lemma is needed to show how to separate n and q 
and afterwards absolutely continuous and discrete part of the Levy measure. 

Lemma 10. Suppose 11 is purely atomic with possible accumulation point only at zero. If f is 
infinitely differentiable away from Gi with locally bounded derivatives and integrable at zero, then 
II* f is infinitely differentiable away from G^+i with locally bounded derivatives. 




Proof. First note that by (|4.15|) for arbitrary a > 
d*" 

(4.19) 

From this simple identity the claim follows for the special case H = Sa- For infinitely many 
atoms the difficulty appears from the fact that 11 is an infinite sum of indicator functions so that 
summation and differentiation in 



(4.20) 



n * f{x) = ( ^ n({a})l[o,a] ) * fix) - J2 n({«}) iMo.a] * f) [x) 



need to be interchanged. To justify the Fubini flip in (4.20) it suffices to show locally uniform 



absolute convergence of the right hand side. An upper bound can be obtained as 

Ui{a}) r \f{y)\dy+ ma}) T \f{y)\dy 

x>aeG x<aeG 

< sup \fiy)\ Y ma})a+ r\fiy)\dy J2 "(W), 

ye[x-c,x] x>aeG "'o x<aeG 

where c denotes the largest atom strictly smaller than x which exists as the atoms are discrete 



away from zero. The right hand side is finite by property (1.1 1, continuity of /, and integrability 
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of / at zero. 

Having justified (4.20), we now show that away from Gi+i 

^'-n * fix) = J2 n({«})i (i[o,a] * /) ix) 



(4.21) 



aeG 



which we prove by showing locally uniform absolute convergence of the series of derivatives. First 
note that as a; ^ Gi+i there is c' > such that {x — c' ,x + c') fl G^+i = 0. As the set B — {a ^ 



G : a > c'} is finite, we appeal to (4.19 1 and the mean value theorem, we obtain for a e S'^ n G 



— ^l[o,a] */(a:) < sup 

"■^ ye[x — c',x\r\G 



-Jiv) « 



Then 



E n(W) 



a£B<=r]G 

< sup 

y^[x—c,x] 



^(l[o,a]*/)(a;) 



dy 



E n({a})a, 



aeB-'nG 



enought to show (4.21) as S is a finite set. 



which again is bounded by property (1.1 1 and local boundedness of derivatives of /. This is 



In total we proved that 11 * / is infinitely differentiable away from G^+i and derivatives may be 
taken termwise. Local boundedness of the derivatives follows from the above estimate. □ 



With the lemmas in mind we can give the proofs of the main results. 

Proof of Theorem^ To prove higher order (non) differentiability, we consider the Laplace inversion 
representation of Proposition j2j for N > — '■ 



(4.22) 



dx^ 



(x) 



d'' 

dx'^ 



1 



dx'^ 2tt 



\dO. 



As N was assumed to be large enough, the integral is everywhere fc-times differentiable by Lemma 
[7] so that the critical points claimed in the theorem have to come from differentiating the finite 
(possibly very large) sum. As motivated before the proof, N can be chosen arbitrarily large as the 
higher order convolutions are smooth enough. 

Since representation ( |4.22 1 is valid replacing k by some n < fc, we proceed by induction showing 
that u^'^^ is fc-times differentiable in x iff a; ^ Gfe. The induction basis for n = 1 is true in complete 
generality (without Assumption (A)) due to Corollary [sj Assume next that the claim is true for 
some n < k and consider (4.22) with fc replaced by n + 1. Since the nth derivative does not 
exist on G„, we only need to consider the (n + l)st derivative on M_|_\G„. The integral term is 
(n + l)-times differentiable as seen above so that we only consider the sum 



i=l 

" (—^Y 

= EW-i*(n+'?r(^) + 



(n+l) 



(x) 



i=n+2 



f l*(n + g)"(x). 



We start with the case q — 0. By Lemma [oj n**(x) is everywhere {i — 2)-times continuously 
differentiable implying everywhere the existence of n + 1 continuous derivatives for the third 
summand. According to Lemma[9]the first sum is everywhere sufficiently differentiable away from 
Gn- Therefore we are left to deal with the middle term. But again according to Lemmajoj n*("+i) 
is n— times differentiable on M4.\G„4.i with jumps of the [n + l)st derivative on G„+i\G,i. This 
shows that the {n + l)st derivative of u exists iff a; ^ Gn+i- In particular, setting fc = n + 1 and 
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using U' 
linearity 



the proof is complete for q — 0. To extend the result to g > 0, we use that by 



n 

(n+,r(.) = 5:(;;)n- 



* q 



■ (n—k) 



(x) 



71— ± 

n*"(a;) + 5] (^) n*^ * q*'-"'^\x). 



fe=0 

Plugging this into the previous equation reduces the problem to the case q = and additional 
convolutions with higher order convolutions of the constant function q. By Lemma [lO] and induc- 
tion, smoothness of higher convolutions of the constant function q implies that the convolutions 
fl*k ^ q*{n-k) j^j.g differentiable away from Gk- As k < n — 1 this shows that the additional con- 
volutions do not generate additional points of non-differentiability. 

Hence, that the claimed differentiability property for u'*' follows from that for u. □ 



Proof of Theorem^ The stra tegy is similar to the one of Theorem [2j Choosing again N > 
— , representation (4.22) holds by Proposition j2j and taking into account Lemma jlj implies 



that we only need to discuss differentiability of H*" = (IIi -I- n2 + q)*" for n — 1, .... N — 1. 
The first term Hi -1-112 -I- is infinitely differentiable precisely for any a; ^ G as Hi is constant away 
from the atoms and jumps downwards on G and II2 is infinitely differentiable. For the iterated 
convolutions we have to separate the contributions of IIi and II2. Writing again 



(Hi + n2 -f qr{x) = J2 Q ^1' * (^2 + qy^^-'H^) 

k=0 



= liTix) + ^ ( J nt'= * (n2 -f 



fe=0 



the problem is reduced to pure convolutions of IIi and mixed convolutions. The pure convolutions 
have been analyzed in the course of the proof of Thcorem[2]and we have the claimed differentiability 
behavior of the theorem. 

The proof is complete if we can show that smoothness of 112 -I- g makes the mixed convo luti ons 
-k) everywhere infinitely differentiable away from Gk 



10 



we 



n^'^ * (n2 -I- *("- *=) e verywhere infinitely differentiable away from Gk. To apply Lemma 
use Lemmajsjand (4.2 1 to see that (112 + is everywhere infinitely differentiable and integrable 
at zero for arbitrary integer I. Hence, Hi * (II2 -I- g)*' is infinitely differentiable away from G and 



furthermore integrable at zero by the same arguments used to derive (|4.2 ) 
shows that H^'^ * (112 



10 



Inductively, Lemma 

' is infinitely differentiable away from Gk for any integers k,l. As 

□ 



k < n — 1 this shows that no additional points of non-differentiability are caused by H2 . 



4.3. Asymptotic Behavior. In contrast to the application to differentiability of the Laplace 
inversion representation, we now apply the series representation to find the asymptotics of u*-'^ 
and it's derivative at zero. The Laplace inversion representation is then applied to the asymptotics 
of u' at infinity. 



Proof of Theorem^ From Equation (3.2) it follows that 



(4.23) 



1 



j,l*{q + n){x) ^ ' ^\^{q + ^V){x 



00 / -I l*(<;-|-n)*"(a;) 



Hence, we need to show that the latter summand of the right hand side converges to zero (abso- 
lutely) as X tends to zero. To this end we use the estimate (4.2) to obtain 



Er=2(-i) 



n+l l*(?+n)-'"(a:) 



5L(i*(g + n)(x)) 



< 



(l*(g + n)(a;))' 
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Clearly, limj-^o 1 * + = by the property ( 1.1 ) so that for x small enough the right hand 

side is bounded from above by 

1 

-, l*{q+n){x) 

^ s 



1. 



This shows that the right hand side vanishes for x tending to zero proving (4.23). The higher 
order asymptotics follow in precisely the same manner. □ 

The refined first order asymptotics can now be deduced directly: 

Proof of Corollary^ By Theorem |4| u'*' behaves asymptotically like ■g2iqx + ^{y) dy) at zero. 
As n decreases we have for x sufficiently small 

u(«)(a;) ^ qx^-x1l{x) _ q + Il{x) 
X ~ b'^x (5^ 
and the right hand side diverges if the Levy measure is infinite. In the finite case we obtain 
similarly for x sufficiently small 



9 + n(a:) ^ ^ g + n(0) 

b"^ ~ X ~ b^ 



proving the claim. 



□ 



Proof of Theorem^ Without loss of generality we assume (5 = 1. Also consider the case when 
n(0+) = oo. To prove the theorem, taking into account Equation (3.2 1 we need to show that for 



(4.24) 



lim J2 i-l)\q + nrH^) = 0. 



k=n+l 



Note that taking into account Lemma [4] we may appeal to q + Il{y) < Cy for e > 

satisfying /3{n) + e < 1. Then for y small enough 



iq + n){z){q + IL{y~z))dz 



By induction, for all y such that q + tl{y) < jg easy to show that for each k there is 

a constant Ck depending only on /3(n) such that 

For £ small enough, the quantity — fc(/3(n) + e) + (fc — 1) is strictly posi tive for fc > n by the 



assumption /3(n) < Therefore lim^_j.o((7 + H) 



and with (4.2) for a; > 



{q + n)*("+l+'^')(x) < sup {{q + n)*("+l)(y))l * {q + Uy^x) 

y<x 

< sup {{q + n)*("+l)(y)) (1 * (g + U){x))' 

y<x 



Again taking into account property ( 1.1 1, we see that lim2:_>.o 1* (q + n)(a;) = and also l im.r^ o(g+ 
n)*("+^'(a;) = 0. Furthermore, we have q + ti{x) ~ ^{x) so that we can easily deduce (4.24). 

If n(0+) < 00 then /3(n) = and clearly (4.241 holds with n = 1. Thus we conclude the 
proof. □ 

For the asymptotic behavior at infinity we need to have a more carefuU look at the Laplace 
inversion representation. 
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Proof of Theorem^ Formally setting A = in the Laplace inversion integral of (3.6 1 leads to 



(4.25) 

(4.26) 
where 



u{x+) ~ 



n(x+) 



^2 

n(x-) 

6^ 



n=2 

N-l , 

y - 

n=2 



1 

2^ 



m" [9) dO, 



D" - 



{-me)) 



1 

" 2^ 

N 



JOx^N 



(5^+1(1 + |^n(0)) 

To see that is well-defined, recall from Remark [s] that for finite mean Levy measure the 
denominator is bounded away from zero. To prove the representation rigorously we need to send 
A to zero in the Laplace inversion representation ( |3.6[ ), ( |3.7| . For this sake, convergence under 
the integral and the interchange of limit and integration need to be justified. Here, the finite 
mean assumption is crucial. As convergence of the exponential is trivial, we only show that 
limA->.o h^{X + id) — {id) = [9) for fixed 9. By definition of it suffices to show 



\\incn*\\ + i9) = m*\ie 



Fir\9) 



for ^ = 1 and I — N . Here we use that the additional assumption E[Xi] < oo is equivalent to 
n(a;) tix < oo to apply dominated convergence justified by the uniform (in A) upper bound 



dx < 



U*\x)dx < U{x)dxj 



< oo. 



Having shown pointwise convergence of h , we now verify dominated convergence for A tending to 
zero. The estimate (4.131 is not strong enough as around the real axis the upper bound explodes 
with A tending to zero. The additional assumption again leads to the (coarse) uniform upper 
bound 



\lm{CU{X + i9)) \ = 



sm{9y)rx{y)dy 



< / n(x) dx < oo 



and similarly for the real part Re{CIl{X + i9)). For large 9 we employ the estimate (4.13) for 
sufficiently small e > so that we obtain from our choice of N the integrable in 9 upper bound 
(uniformly in A) 



(4.27) 



h^{X + i6 



< 



n(a;) dx 



\o\<i, 

\9\>1. 



Exploiting the Fourier inversion representation for N large enough we can complete the proof. 
Apparently, it suffices to show that limj-^oo n*"(x) = Oforn=:l,...,iV— 1 and 



lim 



m^{9) d9 = 0. 



The first is a consequence of I[*"{x)dx < oo for which we take into account (4.2) and the 
assumption. The integral vanishes at infinity by the Riemann-Lebesgue the orem as soon as we 
justify that is integrable in 9. But this follows from the upper bound (4.27) which, being 
uniform in A, is also valid for [9) = h^{i9). 

□ 
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